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Abstract 

We report a test of the universality of free fall (UFF) related to spin-gravity coupling effects by comparing 
the gravity acceleration of the 87 Rb atoms in nip = +1 versus that in mp = -1, where the corresponding spin 
orientations are opposite. A Mach-Zehnder-type atom interferometer is exploited to sequentially measure 
the free fall acceleration of the atoms in these two sublevels, and the resultant Eotvds ratio determined 
by this work is r/s = (-0.2 ± 1.5) x 10~ 5 . The interferometer using atoms in mp = +1 or mp = -1 is 
highly sensitive to magnetic field inhomogeneity, which limits the current experimental precision of our 
UFF test. The work here provides a stepping stone for future higher precision UFF test related to different 
spin orientations on atomic basis. 
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The universality of free fall (UFF) is one of the fundamental hypotheses in the foundation of 
Einstein’s general relativity (GR) [1], which states that all test bodies fall with the same accelera¬ 
tion in the gravitational field regardless of their structure and composition. Traditional verifications 
of the UFF are performed with macroscopic bodies that weight di 
material by torsion balance technique [2-4], free-fall method [5- 
achieved a level of 10“ 13 [3, 


ifferently or comprise of different 
t| or laser ranging mission [8,9], 


4, 9]. There are also lots of work investigating possible violation of 
UFF that may be induced by spin-related forces (see, for example, [10-15]), and UFF tests of this 
kind have been performed with polarized or rotating macroscopic bodies [16- 24]. Here we report 
a spin-orientation related UFF test with quantum objects by atom interferometry. 

JFF tests with quantum objects have earlier been performed with a neutron interferometer 
[25D, and in recent years, were carried out by comparing^the free fall acceleration between differ¬ 


ent atoms or between atoms and macroscopic masses [26[-31]. The motivation of using quantum 


objects is not only for potentially higher precision or associated well defined properties, but also 


for more possibilities to break Einstein equivalence principle on quantum basis [32]. For example, 


the variation of t 


the free fall acceleration with atoms in different hyperfine levels has also been 
tested in Ref. [27] at a level of 10“ 7 . Recently, Tarallo et al. [33] performed an UFF test using the 


bosonic 88 Sr isotope (/ = 0) and the fermionic 87 Sr isotope (/ = 9/2) at a level of 10“ 7 by Bloch 
oscillation. In their experiment, the 87 Sr atoms were in a mixture of different magnetic sublevels, 
resulting in effective sublevel of (m F ) = 0. They also gave an upper limit on the spin-gravity 
coupling by analyzing the broadening caused by possible different free fall accelerations between 
different magnetic sublevels. However, we note that possible anomalous spin-spin couplings (see, 


for example, [20, [34], [35]) or dipole-dipole interaction (see, for example, [36]) between the Sr 


atoms with different magnetic sublevels may disturb, or even cover the spin-gravity coupling ef¬ 
fects in their experiment. Since most models describing spin-gravity coupling imply a dependence 
on the orientation of the spin, we perform a new UFF test with 87 Rb atoms sequentially prepared 
in two opposite spin orientations (Fig. 1), namely m F = +1 versus m F = -1. The correspond¬ 


ing free fall accelerations are compared by atom inteferometry [37-41], which determines the 
spin-orientation related Eotvos ratio [42] as 


ns = 2 


8+ ~g- 


( 1 ) 


8 ++ 8 - 

In Eq.(l), the gravity acceleration of atoms in m F = +1 ( m F = -1) is denoted as g+ (g_) to account 
for possible difference. This provides a direct way to test spin-orientation related UFF on quantum 
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FIG. 1: (color online) Schematic of the spin orientations for 87 Rb atoms in magnetic sublevels m F = +1 
versus m F = -1 of the 5~Si /2 hyperfine levels. The bias magnetic field B defines the external direction 
to which the atoms spin is referenced. And the total angular momentum of each atom (denoted by the F) 
processes around B. 


basis. 


Compared with UFF tests using polarized or rotating macroscopic masses, it is much simpler 
to prepare cold atomic ensemble with pure polarization using stimulated Raman transition [43J]. 
However, with atoms in sublevels m F = +1 or m F = -1, the Zeeman effect is considerable, which 
makes the interferometer highly sensitive to the magnetic field inhomogeneity. What’s worse is 
that the phase shift induced by the inhomogeneity is opposite for atoms in the two sublevels, which 
thus can’t be directly canceled in the final comparison of measured free fall accelerations. For the 
interferometer with atoms in m F , the phase shift induced by the gravitational acceleration and the 
magnetic field gradient is expressed as 


<P\ 


mp 


+k e fig mF + 2 a z ,im F y B (V n + V,-/2 )T 2 


( 2 ) 


where the superscript ± denotes the corresponding direction of k e s in the interfering pro¬ 
cess, with +k eS (— k e s) indicating the same (opposite) directions between k eS and local gravi¬ 
tational acceleration. And T is the separation time between Raman laser pulses, while F eff = 
T yj 1 + 2 t/T + 4 t/tiT + St 2 /ttT 2 is the effective separation time accounting for the effect of finite 
Raman pulses duration ( r is the duration for the n/2 Raman pulse) [37] - In Eq.(2), the second 
term corresponds to that induced by the magnetic field (only the first order of the inhomogeneity 
is considered), where y B is the magnetic field gradient, a z ,i is the strength of first-order Zeeman 
shift for 87 Rb atoms in 5 2 S 1/2 state, V r is the recoil velocity, and V„ is the average vertical velocity 
of the atoms in F = 1 at the moment of the interfering n pulse (in this work, the atoms are initially 
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prepared in F = 1 before the interfering). 

In order to alleviate the influence of the magnetic field inhomogeneity, according to Eq. (2), 
three steps are taken in this work. Firstly, the magnetic field throughout the interfering space is 
mapped [44,[45], and a region where the field is relative homogeneous is selected for the interfering 
to take place, namely making y B as small as possible. The selected region is at about 736 mm 
height above the magnetic-optical trap (MOT) center. And there the magnetic field varies less 
than 0.1 mG over several-millimeters vertical distance, while the magnitude of the bias magnetic 
field is about 115 mG. Secondly, the direction of the effective Raman laser wave number keg can 
be reversed to make a differential measurement for each m F [37]. A majority of the influence 
(the part associated with V n in Eq. (2)) induced by the magnetic gradient will be canceled using 
this differential measurement, since the influence is almost independent off & eff (we note that V n is 
typically much larger than V r ). However, with the Raman lasers configured in +k e ff versus —k e s, the 
directions of the recoil velocities are opposite. This induces a tiny difference between the atoms’ 
trajectories. And consequently causes a residual influence in the differential measurement result, 
which is the part associated with V r in Eq. (2). The third step is to correct this residuum using the 
y B obtained from the common mode result for the two interfering configurations of k e s . According 
to Eq.(2), for each m F , the differential mode measurement result (A if d mp = (Aif Fv - A<p“ F )/2) and 
the common mode measurement result (A ip c mp = (A tf np + A(p~ lF )/2) are respectively 

A K f = -KffgmJl' s - a ZJ m F y B V r T 2 
Aip c mp = 2 a ZJ m F y B V n T 2 

According to Eq.(3), the magnetic field gradient y B can be directly estimated from A(/4 F , as long 
as V n is known. 

In order to maximally suppress the magnetic-field-gradient influence in the differential mea¬ 
surement for each m F , it is required that both y B and V n are the same between the +k e ff and -k e ff 
interfering configurations. The two requirements can be simultaneously satisfied by preparing the 
atomic ensembles in the same average velocity, namely V+ = V~ (V s denotes the average velocity 
of the atomic ensemble after the state preparation, and the superscript + denotes the A: e ff configu¬ 
ration). In this case, the atomic ensembles are in the same region when the interfering is taking 
place, and thus y B is the same, while V n is obviously the same since it is directly determined 
by V s . Though techniques are mature for preparing the atoms state by velocity-sensitive Raman 
transition (VSRT) [43], it is in fact not simple to ensure the equality between V+ and V~. Using 


conventional state preparation method (see, for example, [46]), the equality will strongly depend 
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on the pre-determined Zeeman shift and AC-Stark shift. And the corresponding variations will 
cause opposite changes for V+ and V~. Here we explore an easy but reliable method to guarantee 
this equality. For the two interfering configurations, we implement the state preparations using 
the Raman lasers both configured in +/c e ir with the same effective frequency oj e ir (oj e ir = a>\ - oj 2 , 
namely the frequency difference of the two laser beams in Raman lasers.). In this case, for each 
trip, the state preparations are completely the same for the two interfering configurations, and 
thus the average velocities of the selected atoms V s are naturally the same. Compared with con¬ 
ventional operation of the interferometer, in addition to usual Raman lasers frequency chirp, this 
method needs an extra shift of a> eS after the state preparation. This shift will switch the Raman 
lasers configuration from +k eS to -k eS for the interfering process where the Raman lasers need to 
be configured in -k eS . This can be realized by using two arbitrary function generators (AFG) to 
mix with a microwave signal source in the Raman lasers’ optical phase locking loop (OPLL), with 
one AFG to implement the m ef f shift and the other to implement the chirp. 

The experiment is performed in an atom gravimeter that has been previously reported in detail 
in Ref. [40]. It takes 727 ms to load about 10 8 cold 87 Rb atoms from a dispenser using a typical 
MOT. Then the atoms are launched upward and further cooled to about 7 /iK with a moving 
molasses procedure in the atomic fountain. The apex of the fountain is at 750 mm height above the 
MOT, close to the aimed interfering region at about 736 mm, which is helpful to limit the atoms’ 
flying distance during the interfering process. After a flight time of 324 ms from the launch, a 
Raman n pulse with a duration of 26 /us is switched on to implement the state preparation. With 
the Raman lasers configured in +k e s, the detune (defined here as the difference between oj c h and 
the the hyperfine splitting of the two ground levels (5 2 S 1 / 2 ) of the 87 Rb atom) that selects the 
maximum m F = +1 atoms is found to be -1546 kHz, and that for m F = -1 is -1866 kHz. 
After the unwanted atoms are removed by a blow-away beam, the atomic cloud arrives at 736mm 
height and undergoes the 7r/2 — jt — njl Raman pulses with a pulse separation time of T = 2 ms. 
With larger T, the interferometry fringe would become invisible due to different magnetic field 
inhomogeneity experienced by respective atom in the ensemble. This is the reason why the effect 
of the finite Raman pulses durations must be considered in Eq. (2). The transition probability 
of the atoms after the interfering is obtained through a normalized fluorescence detection when 
the clouds falls back into the detection chamber. The entire process of a single shot measurement 
as described above takes 1.5 s. Before the formal data acquisition, V n should be measured to 
deduce y B from A(p c mp ,. This velocity can be obtained from the spectroscopy of the VSRT [43] 
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FIG. 2: (color online) Fringes for different combinations of m F and k c (f, where each fringe shown is an 
average of 10 fringes with one corresponding combination. In one cycle, the fringes arc obtained in turn for 
the combinations of m F = +1 and +k e ff (black circle), nip = +1 and -k e g- (red square), nip = -1 and +k ct \- 
(olive empty circle), nip = -1 and -k e ff (blue empty square). 

with a Raman n pulse applied at the right moment. Here the spectroscopies with the Raman lasers 
configured in +k e ff and -k c n are combined to make a differential measurement, in which method 
the knowledge of the Zeeman shift or the AC-Stark shift is not needed. The measured average 
velocity is V n = 509.0(1) mm/s for the selected atoms in |F = 1 ,m F = +1), and is V„ = 509.4(1) 
mm/s for |F = 1, m F = -1) at the moment of the interfering n pulse. 

Finally the measurement of the gravity acceleration of the atoms in different magnetic sublevels 
is performed sequentially, with different interfering configurations of the Raman lasers. One full 
interferometry fringe is obtained by scanning the chirp rate of m e ii in 20 steps for each m F in each 
interfering A: e ff configuration, namely 30 s for a full fringe. Meanwhile, in order to reduce the effect 
of possible drift of related quantities, for example, the Raman lasers power, four adjacent fringes 
are grouped as a cycle unit, with one fringe corresponding to one combination of m F and k cif . The 
switches between the combinations are automatically controlled by the computer through tuning 
the Raman lasers detune, and the typical fringes for the four combinations are shown in Fig.2. The 
measurement is repeated about 28 hours from cycle to cycle, and the phase shifts are extracted 
by the cosine fitting from the fringes. The differential mode result and the common mode result 
are obtained from the combinations of the corresponding phase shifts. The Allan deviation for 
the gravity acceleration measurement is calculated from the differential mode result for each m F , 
which is shown in Fig.3, and the statistics of the combined measurement results for each m F is 
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FIG. 3: (color online) Short-term Allan deviations for the gravity acceleration measurements using atoms 
in different m F . The Allan deviations for m F - +1 (black circle) and nip = -1 (blue square) are calculated 
from the differential measurement results, while that for nip = 0 (gray empty circle) is calculated from the 
phase shifts obtained consecutively with Raman lasers always configured in +k cF . 

also acquired as shown in Table I. The Allan deviation for the measurement using the atoms in 
m F = 0 with the Raman lasers configured in +k e ff for T = 2 ms is also shown in Fig.3 as a 
reference. According to the Allan deviations, the short term sensitivity for the interferometers 
using atoms in m F = +1 is about 3.4xl0 _3 g/ VHz, which implies a sensitivity of 2.4xl(T 3 g/VHz 
if only one combination of nip and k e g is consecutively repeated. This induced sensitivity is about 
four times worse than that using atoms in m F = 0, which is most probably due to the fluctuation of 
the location where the atoms interacts with the interfering pulses (and thus the fluctuation of the 
experienced magnetic field gradient by the atoms). This fluctuation is caused by the variation of 
the launch velocity as well as the initial launch position of the atomic cloud from shot to shot. 

The uncertainties in Table I are the corresponding statistical standard deviations. From the 
common mode results, the magnetic field gradients experienced for atoms in each m F are deduced, 
which are nearly equal as expected since the interfering region is the same. According to Eq.(3), 
this gradient with a magnitude of -37/rG/mm corresponds to a +8 mrad residual effect for the 
differential result with atoms in m F = +1, and -8 mrad residual effect for m F = -1. It shows that 
a majority of the phase shift due to the magnetic field inhomogeneity is canceled in the differential 
measurement, and the residual is only about 1.2%. And the residual effect due to the Raman pulses 
durations is far less than that level, which is thus safely neglected in this work. In this differential 
measurement of the gravity acceleration with a rather short separation time, some disturbances, 
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for example, that induced by nearby masses or tilt of the Raman lasers, are common for the atoms 
in m F = +1 and m F = -1 and thus cancel in the final comparison, and other disturbances, for 
example, that induced by the AC-Stark shift or the Coriolis effect, can be neglected at the present 


level of accuracy. The Eotvos ratio is finally given by 

t g+~8- 


n s = 2- 


8 + + 8 - 


y mF =u - vt F =~ i 
vt =+ i 


(4) 


where is the corrected differential results as listed in Table I. The resultant Eotvos ratio 
determined by this work is (-0.2 ± 1.5) x 10 -5 , which indicates that the violation of WEP has not 
been observed at the level of 1.5 x 10 5 for the atoms with different polarization orientations. 

In conclusion, we have tested UFF with atoms in different spin orientations based on a Mach- 
Zehnder-type atom interferometer, and the violation of UFF isn’t observed at the level of 1.5 x 1CT 5 . 
This work represents the first atom interferometer which simultaneously measures the gravity ac¬ 
celeration and magnetic-field gradient, and also presents a direct test of spin-orientation related 
spin-gravity couplings on quantum basis. The present precision is limited by the fluctuation of the 
atomic fountain arising from the atom launch procedure. The influence of this fluctuation can be 
alleviated with a more homogeneous magnetic field in future measurement, and in this situation 
the pulses separation time can also be enlarged, which will effectively improve the interferome¬ 
ter sensitivity. On the other hand, the standing optical waves can be explored to manipulate the 
interfering of the atoms (see Ref. [270, f° r example), in which case the internal state of the atom 
doesn’t change and thus the influence of the magnetic field inhomogeneity is dramatically de¬ 
creased. We anticipate a better result for the UFF test with atoms in different spin orientations 
using interferometers of this kind in future. 


TABLE I: Statistics of the differential mode measurement and the common mode measurement. The 2n 
ambiguity is easily removed in this work thanks to the rather short separation time T. The magnetic field 
gradient is deduced from A ip c mF , and the corrected A <p^ in the last column is the result of subtracting the 
residual magnetic field inhomogeneity effect from the original A ip^ in the second column. 


mp 

A < 

Wn F 

JB 

Corrected A tffj 


rad 

rad 

gG/mm 

rad 

+i 

-644.322(7) 

-0.674(8) 

-37.6(4) 

-644.330(7) 

-i 

-644.339(7) 

+0.666(8) 

-37.2(4) 

-644.331(7) 


8 










We thank Zhifang Xu for enlightening discussions. This work is supported by the National 
Natural Science Foundation of China (Grants No. 41127002, No. 11204094, and No. 11205064) 
and the National Basic Research Program of China (Grant No. 2010CB832806). 


* E-mail: zkhu@mail.hust.edu.cn 

[1] C. Misner, K. Thorne, and J. Wheeler, Gravitation (Freeman, San Francisco, 1973). 

[2] Y. Su, B. R. Heckel, E. G. Adelberger, J. H. Gundlach, M. Harris, G. E. Smith, and H. E. Swanson, 
Phys. Rev. D 50, 3614 (1994). 

[3] J. H. Gundlach, G. L. Smith, E. G. Adleberger, B. R. Heckel, and H. E. Swanson, Phys. Rev. Lett. 78, 
2523 (1997). 

[4] S. Schlamminger, K.-Y. Choi, T. A. Wagner, J. H. Gundlach, and E. G. Adelberger, Phys. Rev. Lett. 
100, 041101 (2008). 

[5] T. M. Niebauer, M. P. McHugh, and J. E. Faller, Phys. Rev. Lett. 59, 609 (1987). 

[6] K. Kuroda and N. Mio, Phys. Rev. Lett. 62, 1941 (1989). 

[7] S. Carusotto, V. Cavasinni, A. Mordacci, F. Perrone, E. Polacco, E. Iacopini, and G. Stefanini, Phys. 
Rev. Lett. 69, 1722 (1992). 

[8] J. O. Dickey, P. L. Bender, J. E. Faller et al., Science 265, 482 (1994). 

[9] J. G. Williams, S. G. Turyshev, and D. H. Boggs, Phys. Rev. Lett. 93, 261101 (2004). 

[10] F. W. Hehl, P. von der Heyde, G. D. Kerlick, and J. M. Nester, Rev. Mod. Phys. 48, 393 (1976). 

[11] A. Peters, Phys. Rev. D 18, 2739 (1978). 

[12] B. Mashhoon, Class. Quantum Grav. 17, 2399 (2000). 

[13] Y. Z. Zhang, J. Luo, and Y. X. Nie, Mod. Phys. Lett. A 16, 789 (2001). 

[14] A. J. Silenko, O. V. Teryaev, Phys. Rev. D 76, 061101(R) (2007). 

[15] W. T. Ni, Rep. Prog. Phys. 73, 056901 (2010). 

[16] Z. K. Hu, Y. Ke, X. B. Deng, Z. B. Zhou, and J. Luo, Chin. Phys. Lett. 29, 080401 (2012). 

[17] H. Hayasaka and S. Takeuchi, Phys. Rev. Lett. 63, 2701 (1989). 

[18] J. E. Faller, W. J. Hollander, P. G. Nelson, and M. P. Mchugh, Phys. Rev. Lett. 64,825 (1990). 

[19] J. M. Nitschke and P. A. Wilmarth, Phys. Rev. Lett. 64, 2115 (1990). 

[20] D. J. Wineland, J. J. Bollinger, D. J. Heinzen, W. M. Itano, and M. G. Raizen, Phys. Rev. Lett. 67, 
1735 (1991). 


9 



[21] L. S. Hou and W. T. Ni, Mod. Phys. Lett. A 16, 763 (2001). 

[22] Z. B. Zhou, J. Luo, Q. Yan, Z. G. Wu, Y. Z. Zhang, and Y. X. Nie, Phys. Rev. D 66, 022002 (2002). 

[23] J. Luo, Y. X. Nie, Y. Z. Zhang, and Z. B. Zhou, Phys. Rev. D 65, 042005 (2002). 

[24] W. T. Ni, Phys. Rev. Lett. 107, 051103 (2011). 

[25] R. Colella, A. W. Overhauser, and S. A. Werner, Phys. Rev. Lett. 34, 1472 (1975). 

[26] A. Peters, K. Y. Chung, and S. Chu, Nature (Lodon) 400, 849 (1999). 

[27] S. Fray, C. A. Diez, T. W. Hansch, and M. Weitz, Phys. Rev. Lett. 93, 240404 (2004). 

[28] S. Merlet, Q. Bodart, N. Malossi, A. Landragin, F. Pereira Dos Santos, O. Gitlein, and L. Timmen, 
Metrologia 47, L9-L11 (2010). 

[29] N. Poli, F.-Y. Wang, M. G. Tarallo, A. Alberti, M. Prevedelli, and G. M. Tino, Phys. Rev. Lett 106, 
038501 (2011). 

[30] A. Bonnin, N. Zahzam, Y. Bidel, and A. Bresson, Phys. Rev. A 88, 043615 (2013). 

[31] D. Schlippert, J. Hart wig, H. Albers et al., Phys. Rev. Lett. 112, 203002 (2014). 

[32] C. Lammerzahl, Class. Quantum Grav. 15, 13 (1998). 

[33] M. G. Tarallo, T. Mazzoni, N. Poli, D. V. Sutyrin, X. Zhang, and G. M. Tino, Phys. Rev. Lett. 113, 
023005 (2014). 

[34] T. C. P. Chui and W. T. Ni, Phys. Rev. Lett. 71, 3247 (1993). 

[35] A. G. Glenday, C. E. Cramer, D. F. Phillips, and R. L. Walsworth, Phys. Rev. Lett. 101, 261801 (2008). 

[36] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. Pfau, Rep. Prog. Phys. 72, 126401 (2009). 

[37] A. Peters, Ph.D. thesis, Stanford University, 1998. 

[38] Ch. J. Borde, Metrologia 39, 435 (2002). 

[39] H. Muller, S. Chiow, S. Herrmann, and S. Chu, Phys. Rev. Lett. 100, 031101 (2008). 

[40] M. K. Zhou, Z. K. Hu, X. C. Duan, B. L. Sun, L. L. Chen, Q. Z. Zhang, and J. Luo, Phys. Rev. A 86, 

043630 (2012). 

[41] Z. K. Hu, B. L. Sun, X. C. Duan, M. K. Zhou, L. L. Chen, S. Zhan, Q. Z. Zhang, and J. Luo, Phys. 
Rev. A 88, 043610 (2013). 

[42] R. Eotvos, V. Pekar, and E. Fekete, Ann. Phys. (Leipzig) 68, 11 (1922). 

[43] K. Moler, D. S. Weiss, M. Kasevich, and S. Chu, Phys. Rev. A 45, 342 (1992). 

[44] M. K. Zhou, Z. K. Hu, X. C. Duan, B. L. Sun, J. B. Zhao, and J. Luo, Phys. Rev. A 82, 061602(R) 

( 2010 ). 

[45] Z. K. Hu, X. C. Duan, M. K. Zhou, B. L. Sun, J. B. Zhao, M. M. Huang, and J. Luo, Phys. Rev. A 84, 


10 



013620 (2011). 

[46] A. Louchet-Chauvet, T. Fa rah. Q. Bodart, A. Clairon, A. Landragin, S. Merlet, and F. P. Dos Santos, 
New J. Phys. 13, 065025 (2011). 


11 



